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ABSTRACT 


This thesis reports on methods for solving traveling 
salesman problems with time-window constraints. Two types of 
time windows are considered: hard time windows, which are 
inviolable, and soft time windows, which are violable ata 
cost. For both cases, we develop several heuristic proce- 
dures, including some that are based on Stewart's [Ref.6 } 


effective heuristics for the traveling salesman problen 


without time-window constraints. In addition, we aevelop 
exact algorithus for each case, which are based on the 
state-space relaxation dynamic programming  metnod OE 
Christofides,  Mingozzi, aad Toth [Ref.5]. Computational 


experience is reported for all the heuristics and algorithms 


we develop. 
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A. OVERVIEW 


Consider a traveling salesman having to visit n cities 
or customers. He starts from a depot and needs to visit each 
of tke other n-1 cities only once and then return to tne 
depot. The ccst of traveling between any pair of cities 
(expressed in terms of distance, time or cost, etc), say 
from city i tc j, is given as Ci Ina cost Matrix C. fhe 
problem is to design a tour through the n cities that would 
Minimize the total cost of the tour. This is known as the 
Traveling Salesman Problem which is a well-known classical 
Operaticns research prcbiem. 

The TSP is called Euclidean when the cities that must be 
visited all lie on the same ¡plane and the cost of travelirg 
between any pair of cities is the Euclidean distance between 
them. 

The TSP is an NP-complete problem [Ref. 1, 2]. All known 
exact solution methods have a trate of growth of the computa- 
tion time which is exponential in np. On tne other hard, 
heuristic solution methods have a rate of growth of the 
computation tine which is a lcw order polynomial in n and 
have been experimentally observed to perrorm well. ror tnis 
reason, there has been an extensive amount of  researcn 
directed at TSP heuristics. 

In this thesis we consider adding time window 
constraints to the TSP. That is, ir t; is the time tnat the 
EN cuan visits city i, then t;must satisfy 1; S tS uz, 
where l; and u, are the specified lower and upper bounds of 
a time window. This problem is not as well studied as the 
unconstrainted TSP, but there have been a few approaches 


used on the probler. 


Psaraftis [Ref. 3] has 
model and 


which are 


[ Ref. 
bcund procedure  wnich is 
Christofides etal. [Ref. 


programming state 
bcunding in£crmation within 
The this 


heuristic algorithms which 


objective or 


tour that visits 


We 


optimal 


window. are given a 


solution procedure for two 


similar to time-window constrained 


Space relaxation 


presented a dynamic programming 
dial-a-ride problems, 


TSPs. Baker 


4^] has presented an exact algorithm using a branch and 


effective for small 


5] 


very ns 


have presented a 
TO 


a branch and bound algorithm. 


dynamic 
procedure compute 


Study is to develop exact and 


will provide an optimal or near 
time 


X,Y 


each city in its given 


depot location, a set of 


co-ordinates for n cities and a set of time windows . 


A common 
problems. A set of 
among several vehicles. 
decomposes into one 


prospective application 


of large TSPS (with or without time constraints), 


our research and computation 


than 30 nodes). 


we consider two 


ccnstraints: hard time 


windows cannot te 


but a penalty cost must be paid when they are. 


can ke defined individually 


differ for early and late arrivals. 


application of the 

customer orders 
Given a partition, 

TSP for each vehicle. 


and in deference to 


violated. 


TSP is in vehicle routing 


must be partitioned 
the probiem tnen 
Because of this 
the difficulty 


we CONE 


to small-scaled problems (iess 
different kinds of tine window 
windows and soft time windows. Hara 


Soft windows can pe violated, 
The penalties 
fcr each customer, and they can 


Generally, the penalty 


for arriving tefore the lower time window bound is much less 


than tne 
Chapter III, 
Chapter IV, 


we present the 


We developed 


TSP anā time-constrained TSP. 


[Rer. OTE e cent 


penalty for arriving 


several Fortran 


heuristics, 


after the upper bound. Im 


hard time window approach and in 


we present the soft time window approach. 


programs for 
For the TSP, 
CCCO and CCADS 


solving the 
we use Stewart's 


For the time 
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constrained ISP problems, we develop some new heuristics, 
scme of which are modificaticn of Stewart's heuristics for 
the unconstrained problen. 4e also developed exact algo- 
rithms for both hard and soft windows uSing Christofides et 
al.'s [Ref. 5] method of state-Sspace-relaxation dynanic 
programming and branch and bound. This is described in 
Chapters IIi and IV. 

Finally, we describe a hybrid of the heuristic and exact 
programs. Ike hybrid uses the overall structure of tne 
exact program, but the upper bounds are obtained with the 


heuristic. This is discussed in Chapters III and IV. 


B. THE TRAVELING SALESMAN PROBLEM 


A tour is a chain which passes through all the a nodes 
and in which the first and the last nodes coincide. A tour 
is also known as a Hamiltonian cycle. 


Let a tour be denoted by t = ne bt ie 2) and 
: n 


the cost of this tour be 


n-1 
C (t) = C. t OQ a 
ge 19 l 
3=1 do A n 1 
Her e ( E A ) is a permutation of the integers 
n 


from 1 toa, giving the order in whica tke cities are 


visited. 
The Traveling Salesman Problen can be defined as 
follows. Given a graph G = { N,A } composed of a set of 


ncdes N, a set of arcs A connecting these nodes, anda cost 
(distance ) Cis associated with each arc (i,j) TnTAds The 
ISP is the protlem of finding the minimum cost tour of tne 
nodes in N. Tre foliowing mathematical formulation of tke 
TSP is from Stewart ;Ref. 6]. 
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MIN > c. SE (1. 1) 
1) 


1j 
2599 
STD 
> A. PUES 1 J = Wee. . eu (1529 
! 1] 
1 
> x = 1 i - dM (15239 
1] s 
J 
- = 1 j 1 (1.4) 
- y= 1 = — n a 
NE 2.54 e 
221 iz 
J#L j* 1i 
Y ZU X M O i 72,92 (1.5) 
ij 1j 
A A 
14) 
X c 20058 for alle ED (1.6) 
1) 
NEO tor all I ¡TA 
1) 
where m "e 
"ONE ae Li arc(i,j) is oa the tour 
13 i 
0 otherwise 
y. are continuous variables that force the final 
1 
: Scluticn to be on the tour 
( i.e. include every node in the same route ) 
u is a number 2 n-1 , and 
n 1s the number of nodes in the set N. 
The constraints (12) anda R) ensure that each node 
will be visited exactly once, while constraints (1.4) and 


(1.5) force the final solution to be a single tour that 
starts and ends at node 1 (depot). Ihis formulation is not 
directiy used in our TSP programs, but is or interest in a 


general discussion ot the problem. 
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C. TSP WITH TIME WINDCW CONSTRAINTS 


The time-cecnstrained Traveling Salesman Problem is a 
variation of the TSP that includes time window constraints 
on the time to visit some of the cities. The hard time- 
constrained TSP is to find the minimum cost tour subject to 
visiting each city within its time window. 

“For the time-constrained TSP model, we define a contin- 
uous nonnegative variable, t; , to be the time that the 
salesman visits city 1. Since the salesman must return to 
City 1 (depot) at the end of the tour, the formulation 
includes an additional variable, t, , the total time 
required to ccmplete the tour. 

We assume that a complete, symmetric, nonnegative 
distance matrix, IC d is krown and that time is a scalar 
transformation cf distance. Thus time and distance may be 


used interchangeably. 


The following mathematical formulation of the TSP with 
time constraints is from Baker (Ref. 4]. 


MIN t - t (1.8) 
n+ 1 
Se t v t 2 c f 1 = 22 dcl il (1299 
3: 1 1i 
[Sto xA e Cae Jc, d.n (1.10) 
i j PN 
DESSINS 3 
a Ot ae c ise, Soe fh. CPS I) 
n+1 i 1i 
l 
T. oS tons UF Mew oo semen: (1.415) 
1 ai 1 
where t = the time that the salesman visits city i 
i 3 
Ixi = the absolute value of x 
c = the shortest time required to travel froz 
1) 


CL yY 1i to City j 
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l = the lower bound on the time window for the 


1 
salesman to visit city i 
by assumption all 1 2 0 
i 
u = the upper bound on the time window for the 
at 


Salesman to visit city i 
üu 2 1 eestor aT 
1 st 

The constraints (1.9) through (1.12) ensure a nonnega- 
tive arrival time at city i, t; , be obtained for eachmema 
(node 2 through node n) on the tour. The constraint (1.9) 
guarantees that t, , the time that the salesman leaves the 
node 1 (depot) wili be the smallest t value. The absolute 
value constraint (1.10) ensures that the arrival times of 
any two city differ by amount of time sufficient to alloy 
the salesman to travel between the two cities. The 
constraint (1.11) guarantees that t,,, the time the saiestuan 
returns to the depot, will be the largest t value. The 
inequalities (1.12) and (1.13) are nonnegativity and the 
time window ccnstraints respectively. 

Unfortunately, Baker's proposed model for the time- 
constrained ISP is very difficult to solve, Lecause 
ccrstraint (1.10) is nonconvex. Therefore, we will not use 
this formulation directly in cur program. 


Figure 1.1 illustrates the nonconvexity of constraints 


(1.10) for cne i,j pair, the example | t= tl EE The 
feasible region for this ccnstraint is the union of two 
disjoint sets. Taken all together, constraints (1.10) 
define 2^ disjoint sets sets where m = (n- 1) (n- 2) 72, 


which are very difficult to work wien. 
We can see that the time-constrained TSP is very 


difrerent from TSP, even in fcrmuiation. 


-—- 
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Example of Nonconvexity of (1.10) 
in Two Dimensions. 
15 


3 


Figure 1. 


A. OVERVIEW 


Many heuristic procedures have been Geveloped ror 
sclving TSP- Our purposes in this Chapter are to examine 
some of the well-known heuristics, to review Stewart's [ Ref. 
6] recent heuristic, and to compare these approximate tech- 
niques cn the basis of efficiency and accuracy on a small 
number of exauples. 

in general, heuristic procedures are categorized by 
three broad classes: tour construction procedures, tour 
improvement procedures, and composite procedures [Ref. 7]. 
Tour construction procedures start with a Single node anid 
successively add nodes till a tour is built. Tour improve- 
ment procedures attempt to find a fetter tour given ai 
initial tour. Composite procedures ccnstruct a starting tour 
from one of the tour construction procedures and then try to 
find a tetter tour using one cr more of the tour improvement 


procedures. 


1. Tour Construction Procedures 


There are many methods available for constructing an 
initial tour.  rrocedures which have been generally used are 


given below. 


a. Nearest Neighbor ( Rosenkrantz et al. (Rer. 3]) 


Step 1. Start with any node as the beginning of 
a- subtour. - 

Step 2. Find the node closest to tne last node 
added to the subtour. Add this node to 


tne current subtour. 
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Step 3. Repeat step 2 until all nodes are contained 
in the tour. Then, join the first and last 


node. 


b. Clarke and Wright Savings ( Clarke and Wright 
[ Ref. 91) 


Step 1. Select any node as tne central depot 
which we dencte as node 1. 


Step 2. Compute savings S. =C o Co, 


COCR = 253... ,nDn.« Lt} 
Step 3. Order the savings from largest to smallest. 
Step 4. Starting with the largest savings on the 
list, subtours are assembled such that the 
next node added has the largest remaining 
savings - provided that a constraint is not 
violated. Once a pair of nodes i and j have 
been linkea, they remained linked . 


Repeat until all nodes have been assigned. 


Here, tne quantity s, is the amount of travel 

1 
saved if node j is visited directly after i, aS opposed to 
having separate trips from the depot to nodes i and j. 


Figure 2.1 demonstrates the procedure for two nodes i and j. 


C. Insertion procedures ( RosSenkrantz et al. [Ref. 


8) ) 


An insertion algorithm constructs a feasible 
tcur by successively adding one node to an existing subtour. 
This procedure takes a subtour of k nodes at iteration k and 
attempts to determine which node not in the subtour should 
join the subtour next (the Selection step). Then it deter- 


mines where in the subtour it should be inserted (insertion 
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Figure 2.1 Concept of the Clarke - Wright 
Savings Heuristic. 


step). Stewart (Ref. 6] presented the following general 


algorithmic structure. 


Step 1. (Initial SurĽtour) 
Obtain a TSP tour for a subset or the nodes 
N'c Nrin G: 

Step 2. (Selection Step) 
Find a node k€N-N‘* tc be added to the 
existing surbtour. 

Step 3. (Insertion Step) 
Choose an arc(i,j) in the subtcur on N°. 
Insert node k between i and j and add 
Kt Or NM 

Step 4. If N = Nt, then stop 
(We have a Hamiltonian cycle). 


Otherwise, return to step 2. 
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There are many variations on this algorithmic 
structure depending on the procedures chosen for executing 
steps 1,2 and 3. 

Wiorkowski and McElvain (Ref. 10], Or (Ref. 11], 
Stewart (Ref. 12] and Norback and Love (Ref. 131] all present 
insertion algcrithms that use the convex hull of the set of 
nodes N for the initial subset N'. Nemhauser and Hardgrave 
(Ref. 14] have shown that there exists an optimal tour for 
every Euciidean TSP in which the relative order or the nodes 
on the boundary of the convex hull is preserved. This means 
that the optimal tour visits nodes on the boundary of the 
cenvex hull in the same order as if the boundary itself were 
fcllowed. 

Further justification for the use or the convex 
hull fcr the initial subtour is shown empiricaily by 
Stewart's (Ref. 6] computational experiment. He compared 
several insertion heuristics both with and without the 
convex hull as the starting solution. The results show that 
all the inserticn aigorithms are improved by the use of the 
convex hull. Some are improved substantially, others only 


noderately. 


Many criteria have been suggested for the selec- 
tion of the ncde to be inserted in an insertion procedure. 


(1) Nearest Neighbor ( Rosenxrantz et al. [Ref. 


Se a o A e r a c— — a e we ee oe See CEP ee ee eee A = o = e 


EN). Choose the node k that 1S nearest a node in tne 
arent tour=- I- e. , find k - argmin c. s.t. jJEN-N*, 
1 
i €N'. ^ j 
cheapest. Tnsertion .( BRosenkrantz et ali. 
meet. 8] ). Choose the node k that may be inserted at 


minimal increased cost. I.€., find 


feeeabgMin C tc -c S. t. n € N-N', l1,j€ N'. 
n im ny 1 
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(3) Farthest Insertion  ( Rosenkrantz et al. 
(kero Poole Choose the node k that is farthest from a node 
current subtour. I.e., find k = argmaxc | Set. j € NNE 
iem y: 


(4) Arbitrary Insertion ( Rosenkrantz et al. 








{ Ref. 


Koa) 


J): Choose node k randomly from among N-N'. 
(5) Ratio Insertion ( Stewart [ Ref. 12] )- 


Choose the node k such that the proportional increase in 


ccst is minimal. I.e., find k = argmin (c tco) / C 
m im nj 33) 


Sete ni € Nem ee) SNe 


(6) Perpendicular Distance ( Wiorkowski and 


Se ae A = oe — a e ee eee O A 


McElvain (Ref. 10]).. Cnoose the node k that is closest to 
an arc in the current subtour. 

(7) katio Times Distance ( Ord Rer. dE 
Choose the node k Such that the product of ratio and 


distance is  mBinimized. I. e. , find 


k » arguin ((< + CAC +) Coe cd 
1 im nj ij in mj 13 


p D € N-N' , LT € NE 


(8) Greatest Angle ( Norback and Love [Est. 
13) ). Choose the node k and arc i,j such that the angle 
formed ty the two arcs (i,k) and (k,j) is a maximum.  I.e., 


find k = argmax angle{ arc(i,0), arc(t.J) ) S.t. D e NENE 
n 
Le) € "Ww". 
The insertion criteria that have been used 
fall into twc categories. [Ref. 6] 


1. Cheapest Insertion 


Insert the node ké N-N' between those two 
connected ncdes i,j € N' that minimize the quantity 


C NS E cc 
ik kj J 
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2. Identical Insertion and Selection 


Do selection and insertion in the same 


step. 


The best known procedures of this type for the TSP 
are the branch exchange heuristics [Ref. 7]. These branch 


exchange heuristics work as follows. 


Step 1. Find an initial tour. This tour may be chosen 
randomly from the set of all possible tours, 
cr it may be generated ky one or the tour 
tuilding procedures above. 

Step 2. Improve the tour using one of the branch 
exchange heuristics. 

Step 3. Continue step 2, until no additional 


improvement can be made. 


For a given k, we define a k-change of a tour as 
consisting of the deletion of k branches ina tour and their 
replacement by k other branches to form a new tour. A tour 
is k-opt if itis not possible to improve the tour viaa 
k-change. In general the larger the value of Kk, the nore 
likely it is that a k-opt solution will be optimal. 
Unfortunately, the number of operaticns necessary to test 


all k exchange 1S proportional to n^ , Where n is the number 


of nodes in the TSP. Due to this complexity, values of k = 2 
and k = 3 are most commonly used [ Ref. 7]. The 2-opt and 
3-opt heuristics were introduced by Lin [hef. 15] and tne 


k-opt procedure, for k23 was presented by Lin and Kernignan 
eet. 161. 

Or [Ref. 11 J has designea a modified  3-opt tnat 
considers only a small percentage of 3-branch exchanges. 
This modified 3-opt called Oropt by Stewart (Rel. om 
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considers only those branch exchanges which are composed orf 
a string of one, two, or three adjacent nodes being inserted 
between two cther nodes in the current tour. By limiting tne 
number of exchanges that are considered in this way,  Oropt 
reguires many fewer calculations than a full 3-opt. 

Stewart (Bef. 6] made an experiment or the convex 
hull, cheapest angle insertion algorithm (CCA) which will be 
discussed in the next section as a stand-alone algorithm and 
with each or the three post- processors. The algorithms are 
designated CCA, CCA2, CCAO, and CCA3 for tne convex huli 
Cheapest insertion stand-alone, with 2-opt, with Oropt and 
with 3-cpt respectively. He drew two conclusions from his 
computational results. First, the 3-opt requires substan- 
tially more time than either the 2-opt or tne Oropt. 
Second, the Z-opt is dominated by the Oropt and tne 3-opt in 
quality of solution. 

In computation time, Oro pt only looks at 
approximately 3n of the n possible 3-opt exchanges on 
each pass. There are n ways to select the first brane 
times 3 ways to select the second branch, and n-2 ays iS 
select the third branch. 

This accounts for the fairly close times for the 
2-opt and Orort. The quality of CCAO solutions dominate 
CCA2 solutions. On tne other hand, there is littie cr no 
difference between tne Oropt and 3-opt in terms or solution 
quality. 

Stewart's main conclusion from the above experiment 
is that the Orcpt performs as well as a 3-opt in a small 
percentage of the computer time required by a 3-opt, and it 
Should be preferred tc both the 2-opt and the 3-opt for 
Euclidean TSP's. 
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3. Composite Procedure 


The basic composite procedure is a combination of 
the tour construction and branch exchange procedures. It is. 
obtained by appending a branch exchange procedure to the 
tour construction algorithm as a post-processor. The proce- 


dure can be stated as follows [ Ref. 17]. 


Step 1. Obtain an initial tour using one of the 
tour construction procedures. 

Step 2. Apply a branch exchange procedure to the 
solution produced by the step 1. 


Stop wnen no further improvement can be made. 


The ccmposite procedure is relatively fast conputa- 
tionally and gives good results [Ref. 18]. 


Be CCAC 
1. Algorithm 


The CCAO algorithm designed by Stewart [ Ref. 60] uses 
the convex hull of the nodes in N for its initial subtour. 
Then it inserts the nodes not currently in the subtour where 
they may be inserted most cheaply (the Cheapest Insertion 
critericn). It selects the ncde k tc be inserted at each 
iteraticn acccraing to how large an angle is formed by the 
two arcs that must be added to the current  surtour 
(Selection criterion) in order to insert k.  Finaily it uses 
an Oropt to make local imprcvement on the tour constructed 
in the first stage. CCAO means Convex Hull, Cheapest 


Inserticn, Angle Selection, Oropt. 


Algorithm : CCAO 


~ 


Input : Number cr nodes, x and y co-ordinates of alil 
nodes. 


Cutput: Ordered iist of tour, tctal cost. 
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Step 1. (Initial Subtoun 
Find the convex hull of the set of nodes N. 
Call the set of nodes on the boundary N'. 
Let the initial subtour be the nodes of N° 
in the same order as they appear on the 
convex hull. 

Step 2 . (Cheapest Insertion) 
For each node me N-N', find 


Lagun tc "c 
m m ay im m 153] 


S.t. i, j E€ N', 21, Jj 2 connecteos 


Step 3 . (Greatest Angle Selection) 
For the next inserticn, select the  ncde 
that  maximizes tne angle between the arcs 
(2 and (mej) over all m € N-N!. 


I.e, find k = argmax anglef (1 ,m),(u; d 
m m m 


S.t. noe NN SS 
Insert k between i and j and add 
Kee tO NT, : : 
step 4 : If N! = N, go to step 5. 
Otherwise return to step 2. 
step 5 ; Apply an Oropt to the current tour. Stop 


when no further improvements can be found. 
End of algorithm CCAO 


2. Example 


Figures 2.2 - 2.4 illustrate the above algorithm on 
the TSP defined as test problem [1] in Appendix A. Ficse 
the convex huil iS generated for an initial starting 
subtour. This subtour consists of nodes 2, 13, 12, 14,5, 15,72 
A solid line marks the boundary of the convex hull in Figure 
292. 
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Pigure 2.2 Initial Subtour and Insertion. 


In step 2, each of the interior nodes (1,3,6,8,9,10,11,16) 
is associated with a pair of connected nodes on the initial 
subtour (the dashed lines in Figure 2.2). In step 3, tke 
dashed lines that form the greatest angle (closest to 1809) 


identify the node to be inserted (node 10 in this example). 


Figure 2.3 shows the problem after the first three 
insertions( ncde 10, node 1and node 8 in that order). 


Notice that Some nodes not in the subtour are associated 
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Figure 2.3 Intermediate Subtour and Insertions. 


with new node fairs. Figure 2.4 shows the final tour for 
stage one. This tour is now passed to an Oropt post- 
processcr. In this case the tour from stage one appears 


from inspection to be optital, and OOropt will find no 


lmprove rent. 


26 





pn———— TL —————Mo— OO ee ee ee i A A + A A O | OSSA Ta 
n e a be ana e, ro a, = gee GD 


` Figure 2.4 Final Tour of CCAO. 





In addition to  CCAO, CCCO (Convex, Cheapest, 
Cheapest, Oropt) «sas been coded for the purpose of compar- 
ison. The only difference between  CCAO and CCCO is that 
CCCO uses the cheapest selection criterion instead of 
greatest angle of CCAO. 

4e used Sedgewick's [ERhef. 19] package wrapping algo- 
rithm for finding the convex hull (initial  subtour). 
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Starting with some point (called the anchor) that is guaran- 
teed to be on the convex hull (Say the one with the suallest 
y co-ordinate), take a horizontal ray in the positive direc- 
tion and sweep it upward until hitting another point. This 
pcint is on the huli. Then start at that point and continue 
sweeping until hitting another point, etc. The package is 
completely wrapped when the first point 1S included again. 
The following algorithm finds the convex hull of an array 
L(1,...,n) of nodes, the node L(n+1) is used as a Sentinel, 
that is, a copy ori the first node which is used to signal 
completion of the procedure. The variable NH is maintained 


as the number of nodes so far included on the hull. 
Algorithm ; Fackage Wrapping 


Input : Number of nodes, x and y co-ordinates or all 


nodes. 


Cutput: Ordered list of convex hull and number of 


nodes included on the convex huil. 


step 1 . (Initialization) 
Find and duplicate anchor. Icas TrA 


NMIN = argmin y s-t i eenn set 
i 


NH = 0, L(n*1) - L(NMIN). 


Step 2 : (Swap nodes NH and NMIN). 
Put last node found into the hull by 
exchanging it with the NHtn node. 
NH = NH + 1. 
TEMP = L(NH). 
L(NH) = L(NMIN). : 
L(NMIN) = TEME. 
Step 3 : (Compute angle) 


Compute the angle from the horizontal made 


28 


by the line between L(NH) and each of the 
nodes not yet included on the huil. 


Step 4 (Find next hull node) 


Find the node whose angle is smallest amcng 
those with angles bigger than tne current 
value of the ‘sweep* angle (the angle from 
horizontai to the line between  L(NH-1) 
and L(NH)). 

Step 4 : Stop when the first point is encountered 
again. I.e., L(n+1) = L(NMIN) + 
Ctherwise, go to step 2. 


End of algorithm Package Wrapping 


We used Sedgewick's Pseudo Angle for finding the 
smailest angle in step 3, which is coded as tne 'THETA' 
function. This function returns a real number between 0.0 to 
4.0 that is nct the angle made by 11 and L2 with the hori- 
zontal but which has the same order properties as the true 
angle. If dx and dy are the delta x and y distances fron 
some node to the anchor node, then the angle needed in this 
algorithm is  arctangent  dy/dx. However, the arctangent 
function is likely to be slow and it leads to at least two 
annoying extra conditions to compute : whether XxX ıs zero, 
and which quadrant the point is in. 

In this algoritnm we only need to be abie tc compare 
angles, not measure then. Thus it makes sense to use a func- 
tion that is ruch easier to compute than the true angle buc 
has the same crdering properties as the true angie. A qood 
pauncerdate for such a function is sinply dy / {ay + dx). 
Testing for exceptional conditions are still necessary, but 


Simpler. 
Function THETA ( Pseudo Angle ) 


Input : dx,dy (delta x and y distances from some 


node to the anchor node). 


Za 


Cutput : Pseudo angle made by L1 and L2 with the 


horizontal line. 


begin 

dx = x(L2) - x(L1) : ax = abs(ax) 

dy Y(L2) TA A abs (ay) : 
^ if ( dx=0 ) ana ( dy=0 ) then t 


070 
dy / (az rra 


else t 
if dx « Q. then t = 20 E 
else if dy < 0 then t = 4.0 + t 
end 


End of function THETA 


Figure 2.2 shows how the hull is discovered in this 
way. We used Sedgewick's Pseudo Angle for finding the 
greatest angle selection point. 

The data for our test problems is given in the 
Appendix. The computational results are summarized in Tabie 
ju As can be seen in Tacle I, CCAO is faster than CCCO on 
the smail-scaled test problems (below 30 nodes ), but CCCO 
is faster than CCAO on the moderately large sized problems 
(over 50 nodes). Generally, the accuracy is almost iden- 
tical in both cases. 

Stewart [Ref. 6] showed that in a larce scaled 
problen, tne CCAO algorithm outperforms any other insertion 
and selection algorithms. Thus, we are highly motivated to 
use a modification of the CCAO algorithm for solving the 


tine-window ccnstraianed TSP. 


TABLE 1 
CCMPUTATIONAL RESULTS OF CCCO, CCAO 


BEEN uc  — |  — | 
| CCCO CCAO | 
| Problem Number Best % CPU 76 CPU | 
Number of Known. Over Time Over, Time | 

i Nodes n Solution Best (sec) Best (sec) | 
| [1] 16 66.6039 0.00 0.0133 0.00 0.0006 | 
7 [2] 22 469.0288 0.00 0.0233 0.00 0.0100 
| [3] Ze 278.4371 0.00 0.0166 0.00 0.06660 | 
| [5] 51 429. 7000 2272 003897 3.94 0.2562 | 
| [6] 46 552.9000 ISLAS 1.54 0.6889 
l | 
l 


| 


* CPU times in seconds on IBM 3033. 


III. THE ISP WITH HARD TIME WINDOW CONSTRAINTS 


>— eee Á 


A. INTRODUCTION 


The first time-constrained TSP we consider is the case 
in which late arrivals are not allowed, and early arrivals 
must wait for the opening of the time window before they can 
begin to service a customer. This is called the hara time 


window case and it is illustrated in Figure 3.1. 
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Figure 3.1 Diagram for Hard Tire Window Case. 


The hard time window case corresponds to military opera- 
tions and to some civilian distribution problems. Meeting a 
deadline is ccnsidered a critical factor in this case. fhe 
Soft time window case will be discussed in the nextChapter. 

Consider a grapk G = {N,A} composed oz a set of nodes N 
and a set of arcs À connecting these nodes. we now define 
some notation to be used throughout cur discussion of the 


time-window-ccnstrained TSP. 


2 


1 = Lower bound on the time window at node i 
E (early allowable arrival time at city i). 
u  - Upper bound on the time window at node i 
3 (latest allowable arrival time at city i). 
d = Time required to spend at node i. 
(service time at city i). 
SPEED = Constant speed at which the vehicle 
travels. 
dist = Distance from i to j. 
1] 

C | = Iravel time frcm i to j- 
1) 

Note : c, - dist, , / SPEED. 

1) 13 
We use gs and c{i,j) interchangeably. 


depot = Depot (home) ncde. 
L= ( L(1) ,L(2) ,---,L(n) )- 
= A tour with n stops visited in the order 
L(1),L(2) ,---,Lín). 
nMTHMME-AAxprE val time at city i. 
WAIT = Waiting time at node 1 for the hard time 


window. 
We also use l(i), u(i), d(i), ARRVT(i), WAIT(i) and 


u d , @ERVD , WAIT interchangeably. 
i i i i 


HEURISTIC SOLUTION TECHNIQUES FOR BARD TIME WINDOWS 
1. Nearest Neighbor 


The following is a Nearest Neighbor heuristic 


Similar to the one used in the unconstrained TSP. At each 


iteration we add a new node tc the end of the subtour. TE 


is the first node that can be visited from the last node of 


33 


the subtour, 


te necessary 


Algorithm : 
Input 


Output 


Step 1 


Step 2 


Step 3 


Step 4 


Step 5 


Step 5 


Step 7 


Step 8 


taking into account any waiting time that might 


due to the lower time window bounds. 


Nearest Neighbor 
: Number of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 


Ordered list of tour, tctal travel time. 


e (Initialization) 
Start at the depot. 
Let i-depot, Nt = {1}. 

- Compute ARRVT for all nodes ké€N-N' if k can 
be visited directly after i: 


ARRVE = max { ARRVE , l} +td+c_ č . 
k 1 l x ik 


STRE ARRAI > ag then stop (' no feasible 


Solution'). 


. If ARRVT X l1 , then cost AS 
k k K k 


Otherwise, COSE ARRI x 


. ( Nearest Neighbor Selection) 
Choose the rode k€ N-N' such that CM is 
a minimum. I.e; find 


k = argmin cost s.t. j €N-N*. 
J 


- (Insertion) 
Insert k after i, add k to subtour N', and 
let 1 = k- 
. If N'-N, go to next step. 
Otnerwise, return to step 2. 
» Compute total travel time, then stop. 
Totai travel time = max { ARE E. ) 


* d * c 
k k,depct 


End of algorithm Nearest Neighbor 
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Tais solution was constructed by starting at the 
depot and moving to the nearest neighboring customer that 
has not yet keen visited as long as the upper bound level 
was not violated. This heuristic may fail to solve the 


problem. 


2. SCCO 


This algorithm is designed for the case when some of 
the nodes do not have time windows. We call these nodes * 
time free ". 

SCCO is similar to the cheapest selection, Cheapest 
insertion metkod for the unccnstrained TSP, except that the 
nedes with time windows are treated dirferently from the 
time free nodes. The nodes with windows are inserted in 
order of increasing upper time window tound. 

The tine free nodes are inserted between these nodes 
by cheapest selection and cheapest insertion, for as long as 
the upper bcund time window will allcw. In the end, a 
Mcdified Oropt is used to improve the solution. 


There is one possible difficuity with this apprcach. 


It may become impossible to reach some or the tine- 
ccnStrained ncdes before their upper  Lound. In this case, 
we must delete some node(s) from the subtour. Whenever we 


See an upper pound that cannot be satisfied, we select a 
node to delete ty the following criteria. 

The first criterion is the width of the time window. 
Hence, time-free nodes are considered first. Then, if 
several nodes in the subtour are tied for the widest time 
window, we select for deletion the node that results in the 
greatest time saved. The algorithm is summarized as 


- 


fcllows. 
Algorithm : Successive Cheapest Cheapest Oropt (SCCO) 
Input : number cf nodes, x and y co-ordinates of all 


a5 


nodes, time windows for all nodes. 


Output : ordered list of tour, total travel time. 


Step 1 . (Initialization) 
Start at the depot. 
Let i=depot, N* = {i}. 


Step 2 . Set k = argmin u, S.t. j€N-N'. 
J 


If k is time free node, then set k - depot. 
Step 3 . Calculate AS 
ARRVT = max ( ARRVT , 1} + dtc. 
k 8 1 i lk 
Step 4 . If REEVE ES uF e then go to step 5. 


Otherwise, select time free node m € N' which 
results in the greatest tine saved ror dele- 
tion. Delete node m from N', go to step 3. 
Step 5 . Add node k to the subtour N'. 
Step 6 . Insert time free node j€ N-N' between nodes 1i 
and k by cheapest insertion and Cheapest 


selection (same as CCCO) until A not 
exceed u . 
k 


If ARRVT < i , then set ARRVT = 1. 
k k l k k 


Step 7 . if N'=N, then go to next step. 
Otherwise, let i - k and go to step 2. 

Step 8 . Apply the Modified Oropt procedure to the 
current tour. Stop when no rurther 1mprove- 


ments can be tcund. 


End of algorithm SCCO 


- 


"Successive" means select the node  successiveiy by 


the smaliest upper bound. in the SCCO aigorithn, if the 
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salesman arrives before the lower bound of the time window, 
adding waiting tine, we set the arrival time equal to the 
lower bound. 

The Mcdified Oropt prccedure for-improving the solu- 
tion is described beiow. This procedure consider only those 
exchanges that would result in a node being inserted between 


two other nodes in the current tour. 
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Figure 3.2 Current Tour before Modified-Oropt. 
Figures 3.2 and 3.3 are helpful to understand how 


the procedure works. In both figures, i,j,k,l, and m are the 


nodes in the current tour. Nodes 1l and m are considered to 
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be adjacent tc node k. A test 1s then conducted to deter- 
mine if node k can be located between two other nodes, such 
as i and j, sc that it results in reduced total travel time. 
into. Cale we Make the appropriate arc exchanges, then 


update tae total cost and route orders. 


10 11 12 
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Figure 3.3 Improved Tcur after Modified-Oropt. 


In this example, the three arcs (1,j), (k,1), anu 
(k,2) are remcved and replaced by (i,x), (j,k), and (1,m). 
When no further exchanges improve the solution, the algo- 


Lithn terminates. 
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Figure 3.4  Subtour in SCCO Procedure. 


Figures 3.4 - 3.6 illustrate the SCCO algorithm for 
the TSP with hard time windows given in Appendix F as in 
test problem [1]. In this problem 10 of 16 nodes have time 
windows. The ctner 6 nodes are time free. 

First, the subtour starts at the depot (node 16) and 
we insert the node with the smallest upper bound (node 12). 
We examine all nodes which could be inserted between 16 and 


12 as long as the upper bound on node 12 is observed. In 
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Figure 3.5 Intermediata Subtour in SCCO Procedure. 


this case there is no such node. Then we select the next 
srallest upper bound (node 14), add it to the tour, look for 
nodes to insert before it, and continue in this manner. NOW 
we have formed the partial tour 16, 12, 95805 11, 6, 173235 
shown in Figure 3.4. 

às shcwn in Figure 3.5, we can insert Paren 
nodes between node 6 and node 3. These insertions are made 


according to the cheapest insertion and cheapest selec 
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Figure 3.6 Final Subtour for SCCO. 


criteria. we do not make any further insertions because 
they would cause a time window violatior at node 3. 

Figure 3.6 shows the final tour for the SCCO 
heuristic. This tour is passed to a Modified Oropt, but in 


this case it will find no improvement. 


- 
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3. SCAO 


This heuristic is identical tc SCCO except for the 


use of the greatest angle selection criterion for the time- 


free nodes, 
Algorithm : 


Input 


Output 


Step 1 


Step 2 


Step 3 


Step 4 


Step 5 
S tep 


Step 7 


instead of cheapest selection. 


Successive Cheapest Angle Oropt (SCAO) 


Number of nodes, x and y co-ordinates of all 
Rhodes, time windows for all nodes. 


Ordered list of tour, total travel tine. 


(Initialization) 
Start at tne depot. 
Let i-depot, N! = {i}. 


Set k = argmin u S.t. j € N- NEG 
J 


If k is time free node, set k = depot. 

Calculate SET 

AREVI = max { ARRVT , 1 ) 400 toc EP 
k 1 a 1 1k 

DM < Hi , then go to step? 


Ctherwise, Seiect time free node m€N'" which 
results in the greatest time saved for dele- 
tion. Delete node m from N', go to step 3. 
Add node k to the subtour N'. 

Insert time free node j€ N-N' between nodes i 
and k by cheapest insertion and greatest 


angle selection (same as CCàÀO) until  ARRVIM 
R 


does not exceed u . 
K 


If ARBVT < 1 , then set ARRVT =i. 
K k k k 

Let i = k. 

If N'-N, then go to next step. 


Otherwise, go to step 2. 
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Step 8 . Apply the Modified  Oropt procedure to the 
current tour. Stop when no further improve- 


ment can be fcund. 
End of algorithm SCAO 


This algorithm is same as SCCO except greatest angle 


selection is used instead of cheapest selection, as in SCCO. 
SACRA 


This heuristic was suggested by Professor Rosenthal. 
It is designed for the case when the widths between the 
upper and iower bounds of the time windows are relatively 
large. 

In this heuristic, the SLACK is the most important 
concept. The SLACK(i) can be defined as the maximum amount 
of time by which arrival at node i can be delayed 
without causing an upper bound to ope violated for a node 
currently on the tour. 

The SIACK function can be defined as a recursive 


function as fcllcws. 


SLACK(L(i)) » min ( u(L(i)) - ARRVT(L(i)) , 
SLACK(L(i+1)) * WAIT(L(i)) 3} 


where 
WAIT (L(i)) - max (0, l(L(i)) - AREVT(L(1)) } 


The first element of this recursive function is the 
difference between tne upper Lound and arrival time at node 
L(i). The second one is the sum of next node's SLACK ana 
waiting time cz node L(i). The minimum of these two 
elements is a possible delay time of the arrival time at 
ncde L(i) withcut violating tne upper bound of ail nodes 
after L (i) in the current tour. 
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The advantage of this recursive function is that it 
is easy to calculate a possible delay time without calcu- 
lating new arrival times £or all nodes after Lí(i). Ihe algo- 


rithm is summarized as follows. 
Algorithm : SLACK 


Input : Number of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 


Output : Ordered list cf tour, tctal travel time. 


Step 1 a (Initialization) 
Start at the depot. Let N' - {depot}. 


Step 2 . Sort the upper time windows. 


U = ( Te at 
S hs UW Succ. ees : 
A n 
Step 3 . Set k - argnin u SS NTC NU 
Jr 


Step 4 . Find a mode L(ISTAR) after wnich node k 
| can Le iaserted in the current sequence, 
if such a node exists. Go to step 7. 
(The criteria by which we determine if aa 
insertion can be mace are given below.) 
Step 5 . If there is no suck place to insert node xk, 
then try to find a node L(ISWAP) in the 
current sequence such that k can replace 
L(ÍSWAP) and L(ISWAP) nas a good chance of 
being reinserted somewnere else. 
Select ISWAP which has the largest time 
Window width among candidates for ISWAP. 
If there is no candidate, then stop. 
( * no feasible solution * ) 
Step 6 <= Do Swap) (add “k EO N', and delete L(ISWAP) 
from N*, and set k = L({ISWAP) ), 
then update slack and arrival times. 


Go to step 3. 


qu 


Step 7 . Select the node which results in the mininmur 
additional travel, i .e, the node k which 
minimizes the following quantity. 

C(L(I),k) * C(k,L(I*t1)) —- c(L(i),L(1*0). 

Step 8 . Insert k after L(ISTAR), and add k to N°, 
and update slack and arrival times. 

Step 9 . If N'=N, then stop. 


Ctherwise, go to step 3. 


End of algorithm SLACK 


This procedure Starts with sorting an array 
Us i Byres = ne into ascending order using a heapsort 
iret. 20]. Ihis u array is used to select a node k in 
ascending order for insertion. Since the upper bound cannot 
be violated, this step is performed. Then find a node L(I) 
after which ncde k can be inserted in the current sequence, 
if a such a ncde exists. 

There are two tests which must be administered to 
determine ir k can be inserted after L(I). First, the 
arrival time at node k if k succeeds L(I), which is called 
TESTA must not be greater than the upper bound u . Second, 
if kx precedes L(it1), taen the resulting delay in arrival at 
Tt), which is called TEST2 , must not greater than 
SIACK(L(I+1)). We can calculate TEST1, TzST2 as foliows. 


TEST1 = Arrival time at node k if k succeed L (1i). 


DA AERVI C) 1(L(1))3 + d(L(233)) + c(L(i) ,k). 


TEST2 = Delay in arrival at L(i+1) if k precedes L(1+1). 


max ([ TESTI, l(k) y) + d(k) + c(k,L(i+1)). 
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If there exists more than one node L(I) after which 
node k can be inserted, we select L{I) according to the 
criterion ot least additional travel time. This additional 


travel time, called TEST3, is given by 
TEST3 » C(L(I),k) * C(K,L(I*l)) * C(IUCDODDUECEOIMEM 


When we insert node k after L(I), we update the 
arrival times and SLACKs. In the updating process, we 
compute updated values of SLACK only for the nodes whose 
SIACK actually changes as a resuit of the insertion. 

If there is no place to insert node k, we call a 
Subroutine  cailed 'TSWAP'. ISWAP tries to find a node 
L(ISWAP) in the current sequence such that k can replace 
L(ISWAP) and L(ISWAp) has a good chance of being reinserted 
somewhere else.  TSWAP uses TEST1 ans IEST2 to find a candi- 
date for ISWAE and then uses a largest time window width to 
select ISWAP. If there exists such a ISWAP, then we do the 
Swap and update SLACKS and arrival times and try to insert 


again. 


C. EXACT SOLUTION TECHNIQUES FOR HARD TIME WINDOWS 


1. State-Space Reiaxation Procedure 


A dynamic programming model of the time-constrained 
TSP has been developed by Christofides et al. [Ref. 5]. we 
applied their approach to compute bounding information 
within a branch and bound algorithm. 

Consider the TS? defined on the graph G = (N,A) with 
the time windcw constraints, where N is a set of all nodes 
Of G, atid A is a set of arcs. Let R(j) be the set o£ nodes 
from which it is possible to go directly to node j. We can 


initially set R(j) » N- {1 I| l; 


(o tad te > Usj, because 


ii 
it is impossible to go directly from nede i to node j:if the 
earlist possitie arrival time at node j exceeds the upper 


time window cf node j. 
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Let f(S,j) be the duration of the least time path 
starting at node 1 passing through every node of SS S'=N-{1} 
and finishing at node j. For a given S and j, we can caicu- 


late a minimumz arrival time in node j as 


rises) = f (S-j,i) +d + c. ]. (3. 1) 
(5, 3) "T" i5 gii] j.i) i ijl 
Then, 
T(S ERIS, j), Lf Ll < IT(5,j) < E 
= 1. SiE T{s,5) < 1. 
J 
with the initialization: 
f((31 4.3) *» c(1. 3) LINE SS Cees 
J 1j j 
= 1, LET EE c 
— ES e E C ; > u? 
13 j 


In equation (3.1) tbe minimum arrival time in node j passing 
through the nodes in the set S can be described as the sum 
of three terms : the first is the duration of the Least time 
path passing through the nodes in the set S-{j} and ending 
in node i, the second is the time required to spend in node 
i, and the third is the travel time frcm node i to node j. 

The £{S,j) can be calculated rcr all subsets S of S' 
and for all nodes j by using equation (3.1) recursively. 
Finally, the optimum solution can be calculated as 

E Rug + T * e J 

Since the computer storage requirements increase 
exponentially witn the size of tne problem, tnis methcd is 
limited to smail pU he total number or £(S,j), when 


S contains k nodes, is k k j since f(S,j) must be calcu- 


lated for all subsets S of St, and since each node in S nust 
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be considered as a possible end-node jj. Therefore the 
storage requirement for £(S,j) in aa node problen, is given 
by [Ref. 21]. 


n-1 7 ee EE ¿272 E 
= ne a o 
k= k ) 


The storage requirements to solve a 22 node probien 
exceed 22,020,096. For relaxing this limitation, 
Christotides et al. [Ref. 5] proposed a state space relax- 
ation procedure which is analogous to Lagrangean relaxation 
[ Ref. 22] in integer programming. The state space associ- 
ated with a given dynamic programming recursion is relaxed 
in such a way that the solution to the relaxed recursion 
provides a bound which could be embedded in a general branch 
and bound method {kef. 23]. We describe Christofides et al! 
method for doing this below. 

Consider the dynamic programming formulation (3.1) 
The state variable in that formulation is (S,j), and the 
Stage is the cardinality of S. Let g(S) be a mapping from 
the domain of (S,j) to some other vector space (g(S),j)- 
Let: 


Hig(S),j) = ( (9g(57-3) ,1) 1 1,€ (S-jN R(j)) ] (3. 3) 


Since we are interested in lower bounds to the TSP 
with time constraints, H(g(S),j) in (3.3) may be replaced by 
any larger set that is easier to compute. Thus,  H(g(5)mmD 


can be defined by tne following equaticn: 
H(g(S),3 = £ (g(S-3),1) | i€E(g(S),J3) ) (3. 4) 
where MHS- Q7R (4) )) = Eglo MDE 


For calculating the lower bound of_the problem, eauation 


(3.1) can be changed to the following equation: 


Tog (S95) yl J) -i za tc j 3 


(g (S= 2) ce aed 59 
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- min 
1€E(9(5),3) 


This gives us: 


f(g(S),j) = 


1 . 


zs co 


With the initialization: 


f (g (j), j) 


| 
a 


Finally, tae cptimum solution 


rin 
1€E(9(N ,1) 


T(g(S),J), Loe 


[Sr(g'),1) *d 
1 


PrGiS sj) 71). + ey Col (9265) 


257 


T(e)43) Su. 
, if Tlg(S),]) 5 1. 
' duos T (g(S),3) 2 uU,- 

J 


if E < EE < n 
rf c ENSIS 
LE m > UR 
1j j 
can be calculated as 
+ Sc e 


The mapping can be selected frcm any separable func- 


tion. Christofides et al. used the following mapping 
munction. 
g (5) = iSl. (3.7) 
Then equation (3.6) becomes : 
II mn d £[ K=1,1.) + do tc,. ] (3:8) 
1e E(k:)) i 1j 
where k= (Sl » 1. 
This gives us: 
ILU if ice su. 
A | Re ise 
J 
= oœ e lf 


With the initialization: 


(1,3) = c(1,j) , 
= E 
B J 
= E ; 


i(k EE. 
J 


PEs oe Ss e < u 
IE E < 1 A 
ric > u? 

13 J 
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Finally, the cptimum solution can be calculated as 


min ftis") ri TEC z 
A UMS i TUE 


2. Additional Condition 


In the previous secticn, we discussed Christofides 
et al.'s state space relaxation procedure which provides a 
lcwer bound on the TSP by reducing a state space in dynamic 
programming. This lower bound is effective in branch and 
bcund only if it is a tight bound. This is Similar to the 
Case in integer programming where the effectiveness of 
Lagrangean relaxation in producing bounds is relative to tne 
integer programming formulation. A redundant state-space 
ccndition can be helpful to get a better bound. For this 
purpose, an additional condition was used by Christofides et 
al. to avoid loops rormed by three ccnsecutive nodes [ Ref. 
5]. This can be done in the following way. 

Let k = "151 Leti (k: 1) ke the duration of the 
least time path fron tae  initiai state to state (K, j) 
without loops formed by three consecutive nodes. Let 
f(k,j,2) be the duration or the second least time path from 
the initial state to state (k,j) without loops formed by 
three consecutive nodes. Let p(k,j,UJ) be the predecessor of 
jon the path corresponding to f(k,j,n). With the atove 
definition, recursion (3.8) becomes: 


T2423, 1) 2. miu 1 £({k=1,1,0) + AMES ND (3.99 
1e E(k,j 1 19 


where o= l 3f p(k-1,320) 53 


2, Ctherwise. 
Inis gives us: 


I(k,j,0, ifl s T(k,j,l) Su. (3. 19] 
l , if TÍk,j,1) < 1, 
J J 


f (k, j. 1) 
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= c SEE TU > u . 
j 


Recursson tor £í(k,)j,2) can be written in the 


fcllowing way. 


Let: 
T(k,j,2) »,min f A A) + ad +c if, (3.11) 
FE B(k,) i 1j 
*p (k,j,1) 
where m = 1, aif p(k-1,1,1)#j 
= 2, ctnerwise. 
This gives us: 
mike i, 2)°= T(k,],2), Pie dks jj, 2) = 0. (3a 12) 
= 1. r if Tk, 3,2) £< l 
= S ; If T(K); > ie 
J 
The initialization is 
EC ld ec (1,1) PLE EN Mc sou. (5913) 
D TJ 1 
= 1l AIL Ga Sad 
i 1i 1 
= œ rec. ^» Uu 
11 1 
and 
(151,2) = > (3.14) 
Finally, the optimum solution can be calculated as 
, Sin fü mTi,1)9* d. tc : (305) 
L e ES uL) X 11 


Since the additional condition can avoid considera- 
tion of a useful lower bound, we considered r(k-1,1,2) in 
recursion (3.9) and (3.11) only when the predecessor of 1 on 
the path corresponding to £(k-1,i,1) is j. If we do not 
consider the second least time path in case of p(k-1,1, 1)=j, 
then f(g(S5),j) dces not guarantee the lower bound of £(S,j)- 

For this example, let's consider a4 node TSP with 
time constraints. Node A is the starting node. Dis the 
tine free node. The lower bcund of ncdeB is 9, the upper 


bound of node B is 11, the lower bourd of node C is 19, ard 
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the upper bound of node C is 21. Suppose service time at 
each node is zero. Figure 3.7 shows an optimal route for 


this prcbiem. 


Oo gees Oe ee OO 





i i i |. | i i 
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X 
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Figure 3.7 Optimal Boute of Four Nodes Problen. 


From equation (3.13) we can get: 


f (1,C,1) 


| 
mE 
vo] 

a 


i 
~J 
. 
e 
~] 


£(1,D, 1) 


~ 


Now applying equation (3.9) recursively with i=1, for k=2 we 


can get: 
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LN2.,B,1) = 21. 


IE ul 219, 


ecg pte) 17.07 


Similarly, for k=3 
£(3,B,1) =a, 
f(3,C,1) =o , 


We can see easily that f(3,D,1) is not a lower bound of 
f((B,C,LD]),D). 


J. ranch and Bound Procedure 


In this section we introduce branch and bound 
enumeration which is used to eliminate subtours in tae solu- 
tion of the state space relaxation procedure. Since the 
state space relaxation procedure 1S a relaxation of tne TSP 
with time constraints, the solution to the state space 
relaxation prccedure provides a lower bound on the optimal 
value of the ISE with time constraints. Any heuristic solu- 
tion can provide an upper bound. We denote some notation to 


explain this algcrithm as follows. 


FLBD = The lower pound, which is tne optimal solut- 
icn to the state space relaxation procedure, 
on the optimal solution to the TSP with 
time constraints given restrictions at the 
current node. 

Z - Current upper bound. 

STACK - Array which represent decision tree. It con- 


tains arc lists which have the same head in 


optimal tour to the state space relaxation 


53 


procedure given restrictions at the current 
node. 
{c' ] = Travel time matrix given restrictions at 

Ui the current node. 

There are two types cf tree search. One is depth- 
first search, the other is breadth-first search (Ref. 24]. 
We used depth-first search since breadth-first search 
required substantially more storage. Depth-first search 
simply means that when a separation is defined, one of the 
nodes created by the Separation is immediately selected to 
be the next  surproblen, and when a ncde is fathomed, the 
enumeration always backtracks to the most recently created 
live node. 

One of the most important requirements of any branch 
and bound algorithm is tight bounds. The closer the bounds 
are to the optimal solution, the fewer nodes must be enumer- 
ated. We used tne SCCO heuristic, which was described in 
section B.2, aS an initial upper bound. The lower bound is 
obtained from equation (3.15). 

To save computing time we need a criterion to decide 
whether or not the branching should be continued. If Fie 
is greater than Z, then the node is fathomed since explicit 
enumeration need not be extended belcw the current node. 
For branching we consider the arcs which have the same head 
node inthe directed graph since each arc must have a 
different head in the TSP solution. If there is no such 
arc, then that solution is a feasible solution. After ail 
nodes of the tree are fathomed, a feasible soiution whica 
has the same value as the upper bound is an exact solution 
to the TSP with time window ccnstraints. 

The fcllowing branch and bound algorithm is used in 


the programs written for exact solution. 
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Algorithm : 


Input 


Eranch and Bound Procedure 


Total travel time of heuristic, travel tine. 


Cutput: Ordered list of tour, total travel time. 


Step 1. 


Step 2 


Step 3 


Step 4 


Step 5 


Step 6 


Cep T 


(Initializaticn) 

Let Z = the optimali solution of SCCO. 

STACK = empty. 

A 

Compute FLBD diven restrictions defined by 
NECE LEDO TZ) go to step 5. 

conato the tree) 

Put all arc(i,j) which have tne same head j 
in directed graph on STACK. 

If there is no such arc, save feasible route 
and update Z = FLBD then go to step 5. 

Let travel time of arc{i,j) which is in the 
top of STACK be infinite, then go to step 2. 
(i.e., C! =œ.) 

de EdCK) 

It STACK = empty, go tc step 7. 

If travel time of arc({i,j) which is in the 
top of STACK is finite, let travel tine of 
da are aa) Pe infinite, then go to 
step 2. (1-e., ct =œ .) 

Otherwise, let NEU time of that arc (i,j) 
be original travel time of that arc (i,j) and 
remove tnat arc(i,j) from top of the STACK. 
Go to step 5. (i-e., c' =c_ ) 

(termination) ES d 

If tnere is a feasible route, then the 
optimal travel time = Z. 


Otherwise, there is no feasible solution. 


End of algorithm Branch and Bound Procedure 
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we present the results of our computational experi- 
ence with the algorithms of this Chapter in Chapter V. 
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IV. THE ISP WITH SOFT TIAE WINDCW CONSTRAINTS 


er ES pom mts 


A. INTRODUCTION 


The second time-constrainted TSP we consider is the case 
in which both late and early arrivals are allowed by paying 
a penalty cost. The penalties are allowed to be different 
for early and late arrivals. The penalty cost is calculated 


as follcws. 


Upper penalty cost = max { 0, upper penalty constant 


X ( arrival time - upper bound ) |]. 


Lower penalty cost = max [ 0, lower penalty constant 


x ( lower bound - arrival tine ) ]- 


In fact, the upper penalty constant is greater than the 
lower penalty constant in most cases. Figure 4.1 may be 


helpful to understand this case. 


This apprcach makes every problem feasible, no matter 
what the time windows are, i.e, even if it is infeasible in 
the hard time window case. This reflects a practical point 
of view, especially when it is possible to save a great deal 
of mileage by allowing a small amount of time window 
violaticn. 

In this Chapter, we considered one unit of cost to be 
the same aS one unit of time. In real world problems, itis 
possible to get a cost by multiplying traveling time by some 
constant. l 

We use the notation lp and i for the lower ana 


upper penalty cost at node k. A 






lower urper--» 
penalty penalty 
<>” ¡COS cost 







l 


| 


A A rr rr rrr Xr 
early lower upper late 
arrival ae bound arrival 

=-=- =-= ==> | time window pen — 


Eor citys 


A A A KA ey ee a A A ip a > A al aas a i Me 
an ——————— A Gee UM" 


i 


Figure 4.1 Diagram for Soft Tire Window Case. 


B. HEURISTIC SOLUTIONS TECHNIQUES FOR SOFT TIME WINDOWS 


1- Nearest Neighbor 


== Se oe a ee 


This heuristic is similar to the hara time windows 
except it takes into account any penalty cost that might be 


necessary. 
Algorithm : Nearest Neighbor 


Input : Nunber of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 


Output s: Ordered list cr tour, total ons 


Step 1 . (Initialization) 
Start at the depot. 
Let i-depot, N*' = {1}, cost = 0. 
di 


Step 2 . Compute AKRVT for all.ncdes ke N-N' 


ARRVI = ARRVWVT td * cC . 
k 1 a 1k 


cost = cost tudo: cu 
k i 1 ik 
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SEDEM ARVI «X le, then cost = cost + lp . 
k k k k k 

If ARRVTIT >u , then ccst = cost + uo . 
k k k k k 


Step 4 . ( Nearest Neighbor Selection ) 
Select the node k € N-N' sucn tnat cost 
is a minimum. I.e. , find f 
k = argmin pores Set. J ENF". 
Step 5 . ( Insertion ) 
Insert k after i, add k to subtour N', and 
iet i =k. 
Step 6 . If N! = N , then go to next step. 
Otherwise, go to step 2. 
Step 7 . Compute total cost, then stop. 
Total cost = ccst * d *.c 
k k k,aepot 
End of algorithm Nearest Neighbor 
This solution was constructed by starting at the 
depot and moving to the nearest neighboring customer that 
has not yet been visited. The term "nearest" is modified in 
the sense that we adda penalty cost to the travel time if 


the time windcw for city i is violated. 
2.  SCCO 


This algorithm is also designed for the case when 
there isa combination of tight tine window nodes and time 
free nodes. The strict observance of the upper bound in the 


hard time windows 1S replaced by a penalty cost. 
Algorithm : SCCO 


Input : Number ot nodes, x and y co-ordinates of ail 
nođes, time windows for all nodes. 


Output : Crdered list of tour, total cost. 


Step 1 . (Initialization) 
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Start at the depot. 
Let i-depot, N*'- (i), cost = 0. 
1 


Step 2. Set X » argnin u. Sst. e NIHU 


If k is time free Bode, then set k - depot. 


Step 3 . Insert node k in the subtour N'. 
Compute ARRVT 
ARRVT = ARRVT +d +c č . 
k 1 D lk 
Step 4 . Insert time free node j€ N-N' between nodes 


i and k by cheapest insertion and cheapest 


selection ( same as CCCO) unti ARRS does 
K 


not exceed u . 
k 
Step 5 . Update cost x ds 
cost = cost d tac un 
k 1 1 1k 


If ARRVT < 1 , then cost » cost * TONE 
K K k K X 


If ARRVT 5» u , ther cost = Cost rua 
k k k k k 


Ke 
N, then go to next step. 


Step 6ra Ict i 
ICZN? 


Otherwise, go to step 2. 


Step 7 . Appiy tae Modified Oropt procedure to the 
current tour. Stop when no further improve- 


ments can be fcund. 
End of algorithm SCCO 


This procedure is also similar to the cheapest 
selection, cheapest inserticn method fer the unconstrained 
TSP, except that the nodes with time windows are treated 


differently frou the time free nodes. The nodes with time 
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windows are inserted in order of increasing upper time 
window bounds. The time free nodes are inserted between 
those nodes ty cheapest selection and cheapest insertion, 
for as long as the upper bound of the time windows will 
allow. 

In the end, a Modified Oropt is used to improve the 
solution. This procedure consider only those exchanges that 
would result in anode being inserted petween two other 


nodes in the current tour. 
3. SCAO 


This algorithm is also designed for the time window 
Set which is composed of some tight time windows and some 


time free nodes. 


Algorithm : SCAO 
input : Number of nodes, x and y co-ordinates of all 
nodes, time windows for all nodes. 


Output : Crdered list of tour, total cost. 


Step 1 . (Initialization) : 
Start at the depot. 
Let i-depot, N' - {1}, cost = 0. 
i 
Step 2. Set k = argmin u. Set: J EN-N'. 
J 
If k is time free node, set k = depot. 
Step 3 . Insert node k in the subtour N'. 
Ccmpute ER 
ARRVT = ARBVT +d +c. 
k 1 1 ik 
Step 4 . Insert time free node j € N-N' betueen nodes 


i and k by cheapest insertion and greatest 
angle (same as CCAO) until aa cess Du 


exceed es 


o1 


Step 5 e Update cost E 
cost = cost t d +c . 
k i 1 1k 


If ARRVT « 1 , then cost cost MEME 
k k k k k 


If ARRVT > u , then cost = cost + up . 
k K k k 


Step 6 . Iet i ka 
If N! 2 N, then go to next stop. 
Ctherwise, go to step 2. | 

Step 7 . Apply the Modified-Oropt procedure to the 
current tour. Stop when no further improve- 


ments can be fcund. 


End of algorithm SCAO 


This algorithm is same as SCCO except a greatest 


angle selection in stead of a cheapest selection in SCCC. 


C. EXACT SOLUTION TECHNIQUES FOR SOFT TIME NINDOES 


In this section we describe a state space relaxation 
procedure, which 1s adapted from Christofides et al. (Rei. 
5), for soft time windows. They only ccnsidered the TSP with 
hard time windows and without time windows. The difterences 
are as follows. Tne waiting cost is replaced by a penalty 
cost to be paid in the early arrival case. Late arrival is 
allowed, but a penalty cost has to be paid. So we have to 
calculate tne duration and the penalty cost on each possible 
path to decide the least cost path in each stage. We denote 
the penalty ccst on each possible path as PC in this 
section. 

Consider the TSP defined on the graph G = {N,A} witn 
scft time window constraints. Let S* be a set of all nodes 


except starting node. Let S ke a subset or S'. Let f (SIM 
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be the cost of the least cost path starting at node 1 
passing through every node of S and finishing at node j. 
Let T(S,j) be the totai duration of a path corresponding to 
ENS. 1). Let p(S,j) be the predecessor of j on the path 
ccrresponding tc f(S,j). Let lp(t) be the early arrival 
penalty cost function and up(t) be the late arrival penalty 
cost function. For a given S and j, total duration of a 


path can be calculated as 


NH 


T (S,j) 59,1). + d. DECOR T]. (4. 1) 


1] 
i. 


where p (S,j) 


In equation (4.1 total duration of the least cost path 
passing through the nodes in the set S and ending in node j 
can be described as the sum of three terms: the first is 
total duration of the least cost path passing  througn tne 
nodes in the set S-{j} and ending in node i, the second is 
the tine required to spend in node i, and tae third is tke 
travel time frcm node i to node j. The dynamic prograamin3 
recursion to aetermine the least cost path may then be 


stated as 


IS tn eieetio-j,1) tod + ¢ + PC j (4. 2) 
d ES jJ 1j 
where T1 -[ I(S-j,i) *d, * c, ]. 
i ij 
PC = 0 icles To <b. 
= Ip (7T) ; if Tl < 12 
- up (T1-u ) ut T4.» u? 
with the initialization: : 
INE ESc , af 1 s.c Su 
1j j 1J J 
= em * Esc.) POET zd ur 
=e oes Ups uu) , X C A > 4? 
13 1j j : 1j J 


Finally, the optimum solution can be calculated as 


min tst 3) ww d * c 3 
PB. i 11 J 
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Since the computer storage requirenents are 
increased exponentially with the size of the problem, this 
method is limited to small problens. For relaxing this 
limitation, a state space relaxation procedure can be used 
same as Chapter Ill. 

Consider the dynamic programming formulation (8.2) 
The state variable in that formulaticn 1s (S,j), and the 
stage is the cardinaiity of S. Let g (S) be a mapping from 
the domain cf (S,j) to sone other vector space (g(S),)). 
Let: 


d(g(S),J) = t (95-3) 223) DOS (4. 3) 


Since we are interested in lower bounds to the TSP 
with time constraints, H(g(S),j) in (4.3) may be replaced by 
any larger set that 1S easier to compute. Thus, “4H (9 (SIN 
can be defined Ey the following equation: 

H(g(S5),3) -» ( (S(3-3),1) 1 1 € E(g(S), 35) } (4. 4) 


where  S-3C€ E(g (S) ,j)- 


For calculating the lower pound of the problem, recursion 
(4.1) can be cnanged to the following equation: 
T(g3(5),93) [ T(g(S-J) o) tun NEL ee (4. 5) 
1 1j 
where  p(g(S),j) = i. 


Recursion (4.2) may be stated as 
EAG(S) 43) nin 54267 AS di tc, + PC] (420 
ij 


(345-3). ije T 
(ESO n *c + PC] (4.7) 


O ij 
where TV = [ TI(g({S-j) 71) "© Gp sce 
| i ij 
E „if 1 <10 EE 
E E 5 MON 


V 


upiri-u ) PIT 
J 


64 


with the initialization: 


f(g(3) ,3 7 C,. moi lA Gus u. 
1j y E j 

= Gips t A LE Ss. SEM 

=C ; HEUD(ÍCHNEUSU.) , if c i > a? 

13 1 J 13 j 


Finally, the cptimum solution can be calculated as 
Din SI. — ]. 
TITEL i i il 
The mapping can be selected frcm any separable func- 
tion. We used a mapping function (3.7), which is proposed by 
Christofides et al., same as Chapter III. Then equation 
(4.5) becomes: 
ISP Ssi= di) + d +c... ]- (4.8) 
JE 1 
where p({S!{,j) = i 


kecursicn (4.7) may be stated as: 


f (ISI, j) =, ninr EAST 1) d +c . + PC 4.9 
: ie Es JOE, La i Sa 
Mere Ti = [ T{jSi-1,i) +d +c ], 
1 1j 
Pe = 0 faites TT < u. 
z= lp -T1)—, 2s Tj P 
= up (T1-u ) lii T1 > a> 
Eers e n alization: J 
EI) =c., aio eC cur 
1j j 1j j 
OS I) , af c. . $1. 
1j j 1j 1j j 
EN E SUPE. UY, 1É C > LB. 
1j Tj Ww 1j j 


Finally, the optimum solution can be calculated as 


tin £ S el +d m ° 
ao NDS D O 
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2. Additional Condi tion 


In the previous secticn, we discussed a state space 
relaxation prccedure which is adapted from Christofides et 
dze| Ref 2059 That procedure provides a lower bound on tne 
ISP with soft tine window constraints. The additional 
conditicn to avcid loops formed by three consecutive nodes 
was used to get a better bound { Ref. 5]. This can be done 
in the following way. 

=: Let k = |S]. Let £(k,j,1) be the cost of tne least 
ccst path £rca the initial state to state (k,j) without 
loops formed Ey three consecutive nodes. Let f (k:): 2) moe 
the cost of the second least cost path from the initial 
State tc state (k,j) without loops formed py three consecu- 
tive nodes. Let p(k,j,m) be the predecessor of j on the 
path corresponding to £{k,j,m). With the above definition, 


equation (4.8) becomes: 


2(Kk,3,0") = [ P(k=1,1,m) 4 d+ coe m UE (4.10) 
1 19 
where sr (kK, gen. ) tok 
ü= Tr A 


2 , otherwise. 


With the initialization: 
T(1,3,1) = C 


13 
and 
AZ c 
RECUEST Che LOE tka.) can be calculated in the rollcwing 
way. Let: š 


T (kJ; 0n) = [T(k-1, 2 m ta € AAC 
a ij 


Ihis gives us: 
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BK, JN" - min DER + dec” t PEE (4.11) 
Ie El 1 


k,j 1j 
where PC = 0 aie. = | ike, lie = Y 
= DENT * (Kk, 32a)), qM UK, 4, n) S^ dd 
J 
= up (T (k,j,2)-u ), if T'(k,j,u) > u, 
J J 
m= 1 , if Pp (k— yao erg 
= 2, ctherwise. 
_ With the initialization: 
Ell, i = Co. ALIS SU (4.12) 
11 1 11 1 
EE NEM uci c S iL 
1i 1 11 11 Jr 
=C t up(c -u ) y Ll GMT u 
11 I1 1i 


Recursion tor f(k,j,2) can be written in the following way: 


f(k: ]j:2) = _ nin PE Ma. E Velie (4. 73) 
DEE 1 1) 
i*p(k, 52) 
nere PC = 0 a JE E E A n). < u 
- lpi(l.-T'(k, jin), if DU (ko jm) < El 


MR ee if T'(k,j,m) > u? 
J J 


EE: D(K-— 123.4) * 


2 , otherwise. 


£3 
il 


With the initialization: 
£(1,1,2) = oo (4.14) 
Finally, the optimum solution can be calculated as 
ESA + an TIC s 4. 15) 
a 25 0 i1 - 
Since the additional condition can avoid considera- 
tion of a useful lower bound, we considered t(k-1,1,2) in 
recursion (4.11) and (4. 13) oniy when the predecessor of i 
ENNUNHe path corresponding to f(k-1,i,1) is j. If we do not 
consider the second least cost path in case of p(k-1,i,1)=j, 


then £(g(S),j) does not guarantee the lower bound of £(S,j) 
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For this example, let's consider 4 node TSP with 
time constraints. Node A is the starting node. Dis the 
time free node. The lower bound of node B is 9, the upper 
tound of node B is 11, the lower bound of node C is 19, and 
the upper bourd of node C is 21. Suppose service time at 
each node is zero, lp(t)j-t, and up(t)-5t. Figure 3.7 shows 
an optimal route for this prcblem.  Frcm equation (5.10) and 


(4. 12) we can get: 
f(1,B,1) » 10, T(17B,1) = 10, 2 IAEA 


f(üi,c,l T9 TUMCM 


14. T4, pil, ee) =e, 


£(1,D,1) * 7.07, T(1,D,1) 5 7-07, IE EDEN 


Ncw applying equation (4.10) and (4.11) recursively with 
1=1, for k=2 we can get: 


I (2,;B8;1) = nmin [ 94.7, 29:84 ] = 29784 
1l &£ {C7 D} 


8 


T (2,8, 1) 


18:14; p(2: B: 0) =D: 


Samriarlvy, 


E(2,C,1) = 19, T(2,C,1) = 14. 10172 Mp( De DD EE 
£(2,D,1) = 17.07, 1(2,D,1) = 17.07 pee 2 ee 
ECE (Ki= 3, 
£43,541) Tite res 98.1 oua 
T(3.8.1) 2 24.14, D(3P DIESES 
Simiiarly, 
£(3,C,1) 2» 39.84, T(3,C,l1) = 24. TÀI 2Mp( oM T 
EOS Dd eo : 


We can see easily that £f(3,D,1) is not a lower bound of 
f ((B,C,D),D). 
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3. Branch and Bound Procedure 








We used the same branch and bcund procedure used to 
elininate subtours in the solution of the state space relax- 
ation procedure in Chapter III.C. 3. The SCCO heuristic, 
which was descrited in section B.2, was used aS an initial 
upper bcund, and the lower bound was obtained rrom equation 
(4. 15). 

We present the results of our computational experi- 
ence with the algorithms of this Chapter in the next 


Chapter. 
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V. COMPUTATIONAL EXPERIENCE 


AD IEST PROBLEMS 


Four sets of test data are used in this thesis. Test 
problem number {1] is taken from Sedgewick [Ref.19: p.3091. 
The other protlems, numbered [2], [3] and [Sar are from 
Appendix 9.1 cf Ellon et al.*'s text [Ref. 21]. These test 
problems are shown in Appendices A,B,C respectively. These 
published prctlems contain node and depot locations, but 
they do not include time windcws. 

We constructed time windows for test problems 
CL11,(2),(3] by first using the CCAO heuristic on the uncon- 
strained TSP. Time windows were then placed about each node 
such that the CCAO route was feasible. The idea for gener- 
ating time windows in this way comes from Baker [Ref. 25], 
who used the unconstrained Nearest Neighbor heuristic as nis 


starting point instead of CCAO. 


The tire window widths were set to varying sizes ranging 
from 3 to 14. Scme of the tine windcws were rairly tight 
while others overlapped. This is in contrast to Bake 
work, where all the time windcws nave width equal 2 units. 

The last prcblem number [4] is the same as test protlez 
[3], except that the time windows were constructed from a 
Nearest Neighbor solution tc the unconstrained traveling 
salesman problem, as in Baker [ Fef. 23: Figure 5598 
displays the CCAO solution for test problem [3] and Figure 
5.2 illustrates the unconstrained Nearest Neighbor solution 
for the test problem [4]. We found a small error in Baker's 
TSP solution for the Nearest Neighbor (Ref. 25], in that ii 
nearest node from node 16 is node 17, not node 13. The 


resulting cost is actualiy higher, it is 312209, noto dE 
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Figure 5.1 Unconstrained Solution Obtained by CCAO. 


Each of the four sets of test data was used to create 
four test problems. The separate instances differed in the 
percentages of time window constraints that were chosen to 
be in effect. The four cases were 100%, 90%, 75%, and 50%. 
We refer to this percentage as the "time-window percentage". 

A candom number generator was used to decide which nodes 
would have tire windows. Test problems for the time window 


Ccnstrained TSP are shown in Appendices G through V. The 
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Figure 5.2 Unconstrained Solution Obtained by 
Nearest Neighror Heuristic. 


fenalty cost factors can be varied depending on real world 
problems. We used 2 and 5 as the lower and upper penalty 
cost factor. Also we set the service time at each node to 0 
to make it easy to construct the time windows. 

The computational results are presented in Tables II and 
IX The figures reported represent results of our test 


runs for each test case. 
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B. COMPUTATIONAL RESULTS 


1. Hard Time Windows 


As nected in section Sy the Nearest Neighbor 
heuristic often cannot solve the problen, because the 
arrival time of the nearest node frequently violates the 
upper bound. However in test problem [4], the results of 
the Nearest Neighbor are the same aS in the unconstrained 
probiem, because this problem itself was constructed by a 
Nearest Neighbor heuristic. 

Generaliy, the SCCO and SCAO heuristic car pe easily 
applied to the hard time window TSP. According to our 
experiments, if the time window width becomes large reiative 


to the travel time between nodes on the optimal uncon- 


strained TSP route, then the lower percentage time window 
problems become more difficult to satisfy. This pnenomenor 
can be seen in test problem [1]. That is because the other 


nodes in the cptimal route for the unconstrained TSP problea 
could be inserted without causing viclation of the upper 
cound. 

The SIACK heuristic takes Slightly more time than 
the other heuristics. It achieved lcwer accuracy in test 
propiems [1] and [3] in the 50 percent time window case. 
The exact algcrithm can find the exact answer in most prob- 
lens, kut when there are fewer windows in erfect, it takes 
mcre computation time. It cannot solve the 50 percent time 


window problems [2] and [4] within 180 seconds. 
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TABLE II 
COMPUTATIONAL RESULTS OF THE HARD TIME WINDOWS 
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BJ SE ew SS ee CPU LO ee OE NES CN eee) ee eee, ee eee Ge eee END 


ame a a RS EE ES EE: 


* CPU times in seconds on IBM 3033. 


2. Soft Time Windows 








All of the methods tested for soft time windows were 
able to find some answer to every problem within reasonable 
computing time, except for two instances with the exact 
algorithm.  Witn the Nearest Neighbor heuristic, the quality 
of solution is not desirable. In general, the lower time 
window percentage problems have lower solution quality. As 
in the hard window case, on test problem [4], the results of 
the Nearest Neighbor heuristic coincide with the uncon- 
strained TSP heuristic, because this problem itself was 
constructed by a Nearest Neighbor. 

As in the hard time window prcblem, SCCO and SCAO 
generally find an optimal solution except for one problem 
with 50 percent time windows. In test problem [1] with 50 
percent time windows, the SCCO and SCAO values were 215.686, 
165.544 respectively. The exact algorithm could not solve 
the two test problems with 50 percent time windows within 
180 seconds. The reason is that the solutions of the state 
Space relaxations have many  subtours and it takes a long 


time to eliminate these subtours. 


With both hard and soft time windows, the results 
are sensitive tc tbe percentage, width and position of tne 
time windows. In most problems, the fewer time windows 


there are, the lower the accuracy of the heuristics. 
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TABLE 
CONPUTATIONAL RESULTS OF THE SOFT TIME WINDOWS 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This thesis has presented some heuristics and exact 


algorithms for the solution cf traveling salesman problen 


with tine windcw constraints. We considered two different 
kinds of time window constraints : hard time windows and 
soft time windows. Hard time windows are inviolable, 


whereas soft windows may be violated at a cost. 

For both hard time windows and soft time windows, we 
developed some new heuristics, SCCC and SCAO, which are 
modifications of Stewart's unconstrained TSP heuristics 
[Ref. 6] CCCO and CCAO. Also for the hard time window only, 
we developed the SLACK heuristic. We also developed an 
exact algorithm ror both hard and soft window using state 
Space relaxation dynamic programming and branch and bound as 
proposed by Christofides et al. [Ref. 5]. 

The procedures were shown to be effective on some moder- 
ately small sized problems. A Nearest Neighbor heuristic 
was also developed, but it was often unable to solve the 
problem with hard time windews, and it found very low 
quality solutions with soft time windows. Tnis experience 
is consistent with the findings of cthers [Ref. 7] who 
determined that the Nearest Neighbor heuristic does not 
perrorm well cn the unconstrained TSP. 

The SCCO and SCAO are generally effective cn most of tne 
Stail sized problems we tested, except for the problems in 
which less than half the nodes have time windows. Further 
research is needed in order to satisfactorily solve these 
problems. Ancther problem difficulty that may require more 
research is dealing with wider time windows. 

The SLACK heuristic which is used only with hard time 


windows is slightly slower than the other heuristics. 


17 


Particularly, in the lower time window percentage problems, 
the accuracy becomes lower. / 

The exact algorithm succeeded in solving 14 of the 16 
test problems to optimality, but it was too slow to use in 
most of the lower time window percentage problems. This 
algorithm's perfcrmance also depends upon the quality of the 
upper bound which is obtained from the heuristic. 
Additional research is needed to reduce computation tine, 
but- a working program for at least some problems has 


resulted from this effort: 
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TEST PROBLEM (1] 


RL ee a amus Se ee Cumüs cun  CURÜD ee a ee eee 


| node X y | node x Y | 
i 1 3 9 { 11 10 13 i 
| 2 11 1 i 12 16 14 | 
i 3 6 8 | Le 15 2 | 
l 3 4 3 | 14 13 16 i 
i 6 8 11 { 16 12 10 i 
| 7 6 4 { i 
i 8 7 4 Í i 
| g 9 7 i | 
i 10 14 5 | 


-— --— E A A A A CUm Gp CARE A A A O ld es a a ls ds A (unb cum cumüD ee CAE ee E cum 


Depot cc-ordinates : (12,10) 


problem source : Sedgewick (Ref. 19] 
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APPENDIX B 


TEST PROBLEM [2] 


ee ee ee A ee ee eee eee eee ee ee 


— a O SS ES a SS ee eee eee 


node x 


2.929 
301 
309 
217 
218 
Z0 
242 
230 
249 
256 
269 


= O Ww O0 NID Uu £5 W Nt = 


anh am 


Y 


217 
258 
260 
274 
267 
267 
249 
2062 
268 
267 
2o 


{ node 


12 
13 
14 
15 
16 
17 
18 
19 
Z0 
21 
22 


X 


25:7 
259 
315 
329 
2318 
329 
267 
215 
303 
208 
326 


M 


242 
265 
233 
252 
2952 
225 
213 
192 
201 
219 
181 


Depot co-ordinates : 


problem source : 


Eilon et al. 
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(326,181) 


[ Ref. 


zu 


APPENDIX C 


TEST PROBLEM [3] 


node X 


Y 


X 


y 


LLJA787 — 77 WS 25 


151 
$59 
130 
128 
163 
146 
161 
142 
163 
148 
128 


= O oO O NW MN ££ W bw — 


md ad 


264 
261 
254 
252 
247 
246 
242 
239 
230 
232 
231 


| 
l 
| 
i 
| 
| 
| 


12 
13 
14 
15 
16 
17 
BG 
19 
20 
Z1 


22 


156 
129 
146 
164 
141 
147 
164 
129 
155 
129 
145 


217 
214 
208 
208 
206 
193 
193 
189 
185 
182 
215 


amn ë Shee AX L0 


Depot co-orāinates : (145,215) 


problem source : 


Eilon et al. [Ref. 
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21}. 


APPENDIX D 
TEST PROBLEM [5] 


ee ee ee ee ee ce ee ee ee ee ee UR eee eee eee ee 


| pode x y | node x y | node x y | node x y | 
| 1 37 52 į} 14 12 42 | 27 OUS NEN UU 5 6 | 
| 2 49 49 j 15 36 16 | 28 43 67 | 41 S310 EM 
| 3 52 64 | 16 52 41 / 29 58 48 ù| 42 21 M 
| 4 20 26 jf 17 7:27 23 | 30 S227 Sei es 5 650 
i 5 40 30 (| 18 147 33 | 31 37 69 I TIN 
| 6 21 47 | 19 13 13 | 32 38$ 46 | 45.30 MM 
| 7 317 63 | 20 57 58 | 33 46 10 q. 46.253205 MI 
|j 8 31.0602 | 21 62 42 ] 34 6151335 "MOI eee 
| 9 52 33 | 22 42 57 | 35 02 63 I 49 TE S eee 
1 10 51 21 | 23 15 57 | 36 63 69 I 49 miom M 
| 11 42 41 | 24 52 | 37 32 22 | "75000500 M 
p r2. exhi 3265 I0 25 1 38 | 383085935 M0 l 
PAs 5 25. | 26 27 68 » 39 5I l 


-— mer ee Áe AA e CD O «O: TUTO A: CD O E CI A A: E ee ee ee ee eee ee ee ee ee ee ee ee CURED Á å ee CURED CURED AGER AA 


Depot co-ordinates : (30,40) 
problem source : Eilon et al. [Ref .21]. 
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Í = amc etos A ee AREE 


| node x y | node x y | node x y | node x y | 
met 22 22 1 20 66 14 | 39 30060 | 58 40 60 | 
| 2 36 26 | 21 44 13 | 40 3050 | 59 70 64 | 
| 3 21 #45 | 22 26 13 | 41 1217 | 60 eu Y | 
mee 4535 31- 23 11 28 | 42 15 14 | 61 SE 6 | 
aM» 55-20 | 24 7 43 | 53 1619 62 30 20 j 
1 6 33 34 | 25 17 64 | 44 2148 | 63 20 30 | 
fee 50 —50 | 26 41 46 | 45 5C30 , 684 15 D -| 
1 8 55 y5 | 27 55 34 | 46 5142 | 65 50 70 | 
5 01767 59 |} 25 35 16 j 47 5615 j 66 57 72 | 
i 10 40 60 | 29 352 26 | 48 48 21 | 67 35 42 |] 
mee 55 65 Gf 20 43 26 | 49 1238 į 68 38 33. | 
NEN :51 15531 31. 76 | 50 31556 | 69 50 ü | 
BEER: 551932 122—553. | 51 29 39 | 70 66 a 
meee oe | So 26 29 1 52 54 38 |{ 71 #59 STE 
EN 1062050001 —34 SO 8074053. 55-57 | 72 35 60 | 
8c 215736 1 35 55 50 | 54 067414 | 73 27 24 | 
| 17 33 44 | 36 54 10 | 55 1070 | 74 40O 20 | 
| 18 SEDES Io. 15. 1. 56 6571 75. 40 37 | 
| 19 62 48 | 38 47 66 | 57 6527 | | 


Depot co-ordinates : (40,40) 
problem source ; Eilon et al.-{ Ref .21]. 
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APPENDIX F 
TEST PROBLEM FOR THE SCCO 


| node X y time window | node x y time window | 


| l(i) uG) | I1(i) a 
i 1 3 3 = = l 11 D d 10 17 | 
I 2 111 1 -= - | 12 16i | 
{ 3 O 8 27 36 | TESTO 2 = = | 
ie 41 25057 MES | 14 13 16 5 13 EM 
Ë 5 5 15 = E ] c9 Z mE = = | 
| 6 8 S 23 | A O - - i 
| 7 6 4 35 43 l i 
| 8 7 4 42 49 | i 
| 9 E 7-58 68 i { 
| 10 14 NES 64 | | 


Depot co-ordinates : (12, 10) 
problem source 
node locations : Sedgewick [Ref. 15] 


time windows : see Chapter V. 
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APPENDIX G 
TEST PROBLEN [1-1] 


a a -— XO eee ee ERO G A a SSD ca ai > ee Se G a l 


| node X y tine window | node x y time window | 
| l(i) u(i) l l(i) u(i) | 
l 1 3 92725 32 į 11 10 13 10 17 i 
l 2 11 1 46 53 | 12 16 14 2 9 i 
| 3 6 LE T 36 { 13 15 2 51 59 | 
I 4 4 DEC q5 p 14 13 16 5 13 
5 Se 15 718 28 | 15 20 12 22 30 l 
l 6 8 11 14 23 | 16 12 10 - - i 
{ 7 6 4 35 43 | | 
i 8 7 4 42 49 | | 
i 9 9 7 58 68 i i 
i 10 14 SE 3 64 l i 


Depot co-ordinates : (12,10) 
Give. 2. 0, CU 2 5.0 
problem source 
aode locations : Sedgewick [Ref. 19] 


tire windows z: see Chapter V. 
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TEST PROBLEM ( 1-2] 


ee ee ee ee ee ce ce ee a RE ee ee ee ee ee ee ee ee ee ee ee ee ee ee AMD ee 


X y time window | node x y time window 


l(i) u(i) | l(i) 


u 


Con HD uu 4 t NH = 


WO 


10 


3 SO 32 
Ta 1 46 53 
6 8 27 36 


11 FUORI 10 
12 16 14 2 
DO s 2 5 


| 

| 

| 
4 3 = = | 14 13 16 5 
S US AIS 28 IS 2 wiz 22 
8 11 14 23 | 16 12 10 = 
6 à 35 43 i 
7 4 42 49 | 
J HA OO 68 


i 
14 o 693 64 i 


17 

9 
39 
13 
30 


Am— woe ë è ë a ë De 


oa Be oe 


Depot co-ordinates : (12,10) 
CL = -2:0 S UF S50 
problem source 
node locations : Sedgewick [Rer. 19] 


tine windows : see Chapter V. 
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APPENDIX I 
TEST PROBLEM [ 1-3] 


ÍT A ue Quo rr ee o A A as 


| node X y time window {| node x y time window | 
| l(i) u(i) | EET (3) ee 
| 1 3 9 = - |] 11 10 13 10 17 l 
| 2 11 1 46 53 | 12 16 14 2 9 | 
i 3 6 8 27 36 Im 13 15 2 51 59 i 
| 4 4 323 7 45 | 14 13 16 5 13 | 
i 5 53 352 18 28 | 15 2 12 22 30 i 
| 6 8 11 14 23 | 16 12 10 = = | 
| 7 6 4 35 43 | 
i 8 7 4 42 49 i | 
| 9 9 7 - - | | 
| 10 14 5 = - | { 


Depot co-ordinates : (12,10) 
CL = 2.0, CU = 5.0 
problem source 
node locations : Sedgewick [Ref. 19] 


tire windows : see Chapter V. 


B 


APPENDIX J 
TEST PROBLEM [ 1-4] 


— ED EP EP EP eee ee ee ee ee ee eee ee eee > m 


On HD WM € UY NO 2 


We) 


X y time window | node x y time window 


i 
| 


l(i) u(i) l l(i) u(i) 
3) Segoe | 31 39] 05 17 
11 1-546 53 | 12 16 14  - - 
6 8 - - NECEM 0 9 59 
i Va MUS E ES ES 13 
5 15 - - | 15 2o - 
8 11 14 23 | 16 12 10  - - 
gs EE | 
7 E | 
9 .7 58 68 i 
i4 5 - - 


10 


Depot Co-ordinates : (12,10) 
CL 22:0 707285710 
Problem Source 
node zocations : Sedgewick (Ref. 19] 


time windows : See Chapter V. 
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TEST PROBLEM (2-1] 


GD A aD ae A A A ee ee ee A Gm A A AA TO A es ee A uA ED A A OD TU AA A E A AT Aum ee cw E A ee ee D mm 


j node X y time window | node x y time window | 
| l(i) u(i) | l(i) u(i) | 
| m 295 772 t25 135 | 12 267 242 170 179 | 
| 2 3913-7259 110 118 INE: 765 193 . 202 | 
| 3773091 2260 102 110 jp 14 315 233 53 67 i 
| 4 217 274 242 250 | 155279 2252 81 89 | 
i 5 218 7278 239 246 | 16 318 252 90 98 | 
| A AA Y | 141 149 | 17 329 224 40 4g 
| 7 242 249 279 286 1 0105920607 213 382 393 | 
i 8 230 262 261 271 | 19 275 192 404 413 
| 9 269 268 206 215 | 20 303 201 432 442 | 
| 10 256 267 200 208 E7052 217 99232:52 i 
EE! 265 257 3 193 22m 326 181 = = | 


A A A A E AL A ME A ce ce ee ee DR A A ONE es ee eee Se ee ee A y 
* 


Depot co-ordinates : (326,181) 
Cs e0 CU "5.0 
problem source 
node locations : Eilon et al. [Ref. 21]. 


tire windows : see Chapter V. 


23 


APPENDIX L 


E A A ms c 


TEST PROBLEM ( 2-2] 


GU Gee O a TU A: aa AAA O O ee ee A A A: O O O ee O O cuum O ee es es ee es es es ee ee es ee O O es es ee ee es es O A O O O es ee ee ee ‘a 


y time window 


1 (i) 


u (i) 


i 


a-— «mb 


— O oO ON DM F&F t N = 


X y time window | node x 
l(i) u(i) | 

295 272 A | 12 267 
SOI 47238 110 118 | 13298259 
309 260 102 110 | 142 955 
217 274 242 250 | 1923215329 
218 278 239 246 |. 169p NE 
282 M267 141 149 | 17 329 
242 249 279 286 | 18599287 
239 202 261 271 { 19159225 
249 268 20685215 | 20885059 
295 267 200 208 | 21599209 
20509257 183 193 | 2259326 


242 
265 
233 
292 
25? 
224 
2513 
192 
201 
217 
181 


170 


2 
81 


179 
67 
89 
98 
49 

393 

413 

442 


i 
| 


PCT" S 5 POL... 


awe —_ 0 


UD DU CURED a eun ee ee umm ee ee ee ees ee ee ee ee es ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee eee ee RUND ee ee GARD ee ee TE unm ee nA 


Depot cc-orainates : (326,181) 
CI *2. 0,200 295970 
problem source 

node locations : Eilon et al. 


time windows : see Chapter V 
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Pret. 


21] 


APPENDIX M 


TEST PROBLEM [2-3] 


B e TU A AED ee ee TU eee ee A A y A a De ee eee 


| node X y time window | node x y time window | 
l(i) u(i) l Poy a(i) 4 
1.729500 272 Was” 135 | 
2 301 25 8 = = | 
3 309 260 - - 
DEM MTISEE 71 22 250 
s ms E 246 | 
6 282 267 141 149 E 55329 224 40 49 
7 242 249 - - | 
e]. cene) PAS = = i 
9 249 268 206 215 | 
Q10 256 267 200 208 
112—265. 257 183 193 | 22 326 181 - ~ | 


Depot co-ordinates : (326,181) 
CL = 2.0, CU = 5.0 
problem source 
node locations : Eilon et al. [Ref. 21]. 


tire windows : see Chapter V. 
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TEST PROBLEM ( 2-4] 


| node X y time window | node x y time window j| 
i E(u) | 1(1) UA 
i 1—:295 21272 = - | 12 267 242 170 IE i 
| 2! 230101952598 110 118 | 0713 982 5D = = i 
| 3 309 200 = — |] 14 315 233 57 67 ] 
| 4 217 274 242 250 | 15055: 002 = = | 
i o 279 273 239 246 | 16 7318252 90 98 | 
| 6 282 267 141 149 | 17 329 224 = = | 
i 7 242 249 2797286 | 19952606705 213 = = j 
| Bog) T262 = = | 19 275 192 COI i 
| 39 249 268 206 215 | 20 303 201 432 9000 NN 
{ 10 -256 2267 = = | 21520998217 = = i 
har HA 265. 257 = = | 229832098051 - = | 
Depot cc-ordinates : (320,181) 


CLomrL2.0,-0B = 520 
problem source 
node locations : Eilon et apre cia 


tire windows : see Chapter V. 


eo 


APPENDIX O 
TEST PROBLEM [3-1] 


i| node X y time window { node x y tine window | 
| l(i) u(i) | (e) ue) el 


159 261 185 193 | 13 129 214 259 271 | 


| 

CoS 222 234 | 

163 247 175 185 | 
246 | 1402 154 | 17 147 193 54 68 | 

161 242 166 173 | 

142 239 131 142 | 

163 236 159 165 | 
AS ¢ 21 #139 #182 40 53 | 


=- O0Ó0OoO A WV fF uN = 
ond 
d= 
o 


hb om 


Depot co-ordinates : (145,215) 
Cm = 2.0, CU S 5.0 
problem source | 
ncde locations : Eilon et al. {Ref. 21]. 


time windows : see Chapter V. 
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APPENDIX P 
TEST PROBLEM [3-2] 


| node x y tine window | node x y time window | 
| Leay eae) | 1(1) uN 


E A O TR A D RED AO AO A O A AM O O A A A O A O A A O CHEER O O E E: A A O E O O A AO E A A O O e A a a IT O A A O O O CADO 


| i 

159 261 185 193 | 
130: 254 22919200] 050 
128. 252 T222 2508 i] | 
163 207 174 3185 | i 
240 142 154 | 17 147 193 54 68 | 

161 242 166 173 | i 
142 239 131 142 | | 
| | 

| | 

l 


148 232 1253088 131 


=- O o ON OD 0 F&F t N = 
od 
4-— 
¡en 


ed = 


Depot co-ordinates : (145,215) 
EL SUL 0 -€U S— se 0 
problem source 
node locations : Eilon et al [Ref. 21). 


tire windows : see Chapter V. 
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APPENDIX Q 


TEST PROBLEM [3-3] 


151 
159 
130 
128 
163 
146 
161 
142 
163 
148 
128 


ma O OO OND OO uU 4 WwW N 2 


"ò ò 


264 
261 
254 
252 
247 
246 
242 
239 

56 
232 
Zn 


time window 


l(i) u 


196 


(i) 


204 


y tiwe window 
1 (1) 


217 
214 
2 08 
208 
206 
193 
193 
189 
185 
182 
23 


105 
2 
22 
10 
54 
79 
67 
YU 


u(i) 


116 
10 
105 
19 
68 
89 


Depot co-ordinates 


CL 


. CU - 


problem source 


520 


node locations : 


tire 


windows 


| node x 
i 

p 12 156 
iS ze 
i 14 146 
i 15 164 
i 16 141 
i 17 147 
i 18 164 
ip pice 129 
p 20:155 
EZ — 139 
| 22 145 
(145,215) 


Eilon et al. 


see Chapter V. 
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( Rer. 


20 


APPENDIX R 
TEST PROBLEM [ 3-4] 


| node x y tame window | node x y tiue window | 
i l(i) u(i) | l(i) wie 
| 1 151 264 = - | 12 156 217 105 TEPIS i 
| 2. $159.26 LS TONES | 13 129 214 = = { 
i 3 130 254 = = j 14 146 208 Z 10 | 
| 4 128 252 222 234 | 15997600 - a i 
| 9 163 207 174 185 p 16 141 206 10 19 | 
l 6 146 246 142 154 | 17 147 193 = = | 
i 7 161 242 166 173 | 18 164 193 = = | 
| 9 1142 2839 = = | 19 129980189 30 38 | 
| vf FIGS 236 159 165 | 20 413595 67 75 | 
100 PS: 232 = = { (21 I O = = | 
09 20390532987 2731 = = | .22 09045205 = = | 


O E: A O ETA a EE EE EP a A Ea ES Se SE SS SS SP SS eS ee ee ee ee 4m ONNEND ee m 


Depot co-ordinates : (145,215) 
CLAS O UA 50 
problem source 
node locations : Eilon et alii ker = i: 


tire windows : see Chapter V. 
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APPENDIX S 
TEST PROBLEM [4-1] 


| node X y time window | node x y time window | 
| l(i) u(i) | O uta) 1] 
De: EST 171 179 DEEP EUECM217 72 9 | 
| .2 159 261 162 170 | 13 129 214 237 245 | 
| 3 130 254 196 203 | 14 146 208 5 9 
E 4 1287 252 — 198 206 1 15 164 208 61 #467 | 
| 5 163 2407 128 136 | 16 141 206 10 — 14 i 
6 146 246 106 113 | 17 147 193 22 28 I 
I 7 161 242 122 130 | 18 164 193 48 53 i 
| 8 142 239 97 105 | 19 129 189 261 269 | 
Mero isa 230 138 146 | 20. 155 185 35 40 | 
MO 148 232 89 96 į 21 139 182 273 280 | 
Wii 128 231 220 227 | 2227 145 215  - - 


SU A O CU O ee ee ee UND ee ee ee ee es ee UNDO CURED ee ee ee ee ee ee es ee es ee ee MED ee ee 


Depot co-ordinates : (145,215) 
Cr 220. CU = 5.0 
problem source 
node locations : Eilon et al (Ref. 21]. 


time winaows : see Chapter V. 
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APPENDIX T 
TEST PROBLEM [4-2] 


| node X y time window | node x y time window | 
| l(i) u(i) | l (i) u(i) | 


159- 264 162 170 | 13 129 214 = = | 
130 254 19/0 82 95 | | 
125 27292 196 206 | 15 164 208 61 67 i 
163 247 128 136 | 16 141 200 10 14 | 
246 100 113 i 17 AROS 22 2 | 
| 18 Tiom; 48 53 | 
| 19 129 189 260 Imm i | 
| 20 5 5X OS 35 40 | 
|. 21 Piso N52 ES = | 
128 2351 220.7227 | 2250055125215 = = | 


= OQ uo Oc -J Oo Uu 4&4 WH N = 
sb 
T 
e 


«mb uh 


—] O A a TS A > A TU A A TU A A TU ee A A AAA A O > A A A O TU O A A A O O A A O A O TU A > TO A A TU A A TU O O O A A e 


Depot co-ordinates : (145,215) 
CLAS 2202990225199 
problem source 
node locations : Eilon et al. [ Ref. 21]. 


time windows : see Chapter V. 
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APPENDIX U 
TEST PROBLEM [4-3] 


] node X y time window | node x y time window |] 
i l(i) u(i) | l(i) u(i) | 
12571565 217 72 79 
13 129 214 - - 
130 254 - - 14 146 208 5 9 


151 264 171 179 l 
i 
] 
PAZO Z 198 206 15 164 208 61 67 | 
l 
| 
| 
l 


159 261 = = 


1 | 

2 | 

3 l 

~ | 

5 163 247 128 136 | 16 141 206 10 14 

6 146 246 106 113 i 17 147 193 22 28 

7 161 242 - - |] 18 164 193 48 53 

8 142 239 - - | 19 129; 189 - - 

9 163 236 138 146 O 155 185 35 | 40 | 
0 148 232 89 96 | 21 139 182 273 280 | 
1 


l 
1 
| 
! 
l 
| 
í 
' 
' 
| 
| D 5:099 »2)927 —19922. 145 215  - 2 i 


Depot cc-ordinates : (145,215) 
C= 2-0, CU = 5.0 
problem source 
node locations : Eilon et al. [Ref. 21] 


time windows : see Chapter V. 
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APPENDIX V 
TEST PROBLEM [4-4] 


-—————— = ål å [D ee Ap ee ee e A O a A A A a ee ce De ee ee ee ee ee ee ee ee ee ee 


| node x y tine window | node x y time window | 
| l(i) u(i) | l(i) 1 (IN 
n a diir aa e n | 12 156 217 72 WS 
| 2 159 261 152 170 II GC os a 
DS IO sh = = | 31000598205 5 o 
| 4 128 252 198 206 I TPE O 2 
| 5 163 247 128 136 | 46 12m ON 
| 6 186 246 106 113 | 12 SNR C | 
| 7 161 242 122 30 ā i p (<S O 2 l 
Pe o SIS = z | 19 129 189 261 2599 
| 9 163 236 138 186 | 20 155 185 55 5 D 
Tu A. - E ( 21259980 5 NS - 
eth eto ee 2 a | 22. 46S ee -- - i 


Depot co-ordinates : (145,215) 
I um. UMS erg 
probiem source 
node locations : Eilon et al. [Ref. 21]. 


time windows : see Chapter V. 
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